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Goldstone's theorem. Elimination of the auxiliary field sector at the level of the functional provides 
with 0{N'') accuracy the renormalized effective action of the model in terms of the original variables. 
Alternative elimination of the pion and sigma propagators provides the renormalized NLO effective 
potential for the expectation values of the A^- vector and of the auxiliary field with the same accuracy. 
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I. INTRODUCTION 



<~i ■ Large- TV expansion is a classical non-perturbative tool of quantum field theory [U, 0, Q . Leading order solution of the 
Oh' 0{N) symmetric model has been applied to interesting problems of finite temperature phenomena, in particular the 
study of the restoration of the SU{2) x SU{2) ^ 0(4) chiral symmetry of QCD [1, d, Q. One of its attractive points 
is that it preserves Goldstone's theorem at every order of the expansion, a feature not shared by all resummations 
of the original perturbative series. Being a resummation to all orders, it has some features which are absent at any 
' given order of the perturbation theory, but which are believed to be true for the exact solution. One example is the 
I renormalization scale invariance observed at leading order (LO) of the 1/A^ expansion. On the other hand, related to 
J> ' the now well established triviality of scalar theories (see e.g. for a review), it shows the presence of a tachyonic pole 

[ (Landau ghost) in the renornialized propagators. 
C I In early publications 0, [sl, |^ the appearance of tachyons was considered as an inconsistency of the large- TV approx- 

■ imation. In next-to-leading order (NLO) investigations of the expansion, started already in 1974 Q, the tachyonic 
, problem seemed to be aggravated because the renormalized effective potential appeared to be complex for all values 
' of the field. This led to the claim that the large- TV expansion breaks down. Extensive studies of the 0{N) symmetric 
I model [1, [13] established the consistency of the 1 /TV expansion for the effective potential and revealed its rich phase 

■ structure. Considered in a restricted sense, as a renormalized effective theory, the large- TV expansion turned out to be 
' a valuable nonperturbative tool, when applied to phenomena dominated by scales much lower than the cut-off. A strict 

J> . cut-off version of the model was considered in [ll| showing that with restrictions on the value of the background field 
the model has a phase with spontaneously broken symmetry, free of tachyons. A different solution to the tachyonic 
problem, called the tachyonic regularization, was proposed in [l^.[l3|. where the tachyonic pole is minimally subtracted. 
Most of the studies were realized in a reformulation of the model in which the quartic self-coupling of the TV-component 
scalar field is replaced by an auxiliary field mediated interaction. 

The renormalization of the 0{N) model was performed at NLO in Q. While pointing out all the divergent integrals, 
this analysis skipped the explicit calculation of the counterterms which were determined only at LO. Calculation of 
the self-coupling counterterm in [TJ] showed that at NLO the /3- function is corrected by terms of order 1/TV. The 
success of the renormalization program of the two-particle irreducible (2PI) approximation (Tsl . ITgI . [iTj revived also the 
study of the renormalizability of the NLO corrections of the large- TV expansion [l^, [l^. Despite all these efforts, the 
detailed and explicit knowledge of the counterterms is missing in the literature, even today. More recent publications 
with interesting finite temperature applications raised doubts about the renormalizability of the NLO approximation 
for arbitrary values of the field expectation and away from the saddle point value of the auxiliary field [lo, [2l| . This 
issue was addressed in an original paper by Jakovac |22l | , where a momentum dependent counterterm is introduced for 
the auxiliary field, which keeps its "propagator" at its classical expression. This generalization of the set of allowed 
counterterms was shown to lead to renormalizability for arbitrary field and auxiliary field expectations. 

The present contribution provides an explicit construction with 0(TV'^) accuracy of the generalized renormalized 
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effective action in the auxihary field formulation of the 0{N) model. This generalisation considers beyond the fields 
also their 2-point functions as dynamical variables of the action functional. This feature is common with the approach of 
2PI derivable) formalism. Subsequent eliminations of selected partial sets of field variables and two-point functions 
lead to functionals depending on remaining variables. We start from the IPI Dyson-Schwinger (DS) equations of the 
model in the auxiliary field formulation. An action functional will be given for that specific closure of the infinite set 
of DS-equations which keeps the coupling of the auxiliary field at its classical value. This closure condition is called 
sometimes the Bare Vertex Approximation (BVA truncation) p3j . 

The leading order part of the effective action is proportional to the number of degrees of freedom N. At next-to- 
leading order (NLO) our goal is to construct the renormalised effective action with C(7V°) accuracy in an arbitrary 
field background. This goal will be achieved by expanding the propagators in inverse powers of N and omitting from 
the action functional all terms which contribute beyond NLO. Then the actual task is to construct the renormalizing 
counter functional for this approximation, which preserves also Goldstone's theorem, obviously obeyed in the NLO 
expansion of the bare theory. We remark here that the 2PI-1/N approximation where one looks for a self-consistent 
solution of the stationarity conditions of the effective action with respect to the propagators represents a different 
resummation strategy. This is reflected also by the fact that the counter-functional we determine is not applicable in 
the 2PI renormalisation program. 

The fact that 2-point functions will be determined following the hierarchy and not self-consistently has the 
definite backward effect of implying secular behavior in time dependent applications 23, 24]. ft is, however, essential 
for ensuring Goldstone's theorem, since an 0{1/N) expansion on the level of the 2PI generating functional was shown 
to not cure the violation of Goldstone's theorem by the self-consistent propagator [l^. Proposition for the preservation 
of Goldstone's theorem for the self-consistent (2PI) propagators exists at present only at LO (Hartree-level) [1^ and 
follows efforts initiated in the framework of non-relativistic many-body theory [26^ . 

By explicit construction of the counter functional in the auxiliary field formulation we demonstrate that the model 
is renormalizable at NLO in the large- expansion for arbitrary vacuum expectation values of sigma and auxiliary 
fields. Elimination of the auxiliary field and related propagators at the level of the 0{N'^) accurate functional leads to 
the recovery of 0{N) and 0{N°) terms of the NLO 2PI effective action of the 0{N) model [27], this time completed 
with all renormalizing counterterms of the same accuracy (valid, however, only in case of NLO large N expansion of 
its propagators). Alternatively, elimination of the NLO pion and LO sigma propagators produces the renormalized 
effective potential for the sigma and the auxiliary fields i2Qi], with correct counterterms. 

The method of the analysis and separation of the overall and subdivergences follows very closely our previous work 
[2M , on the counterterm construction to 2Pl-functionals in constant field background. Here, the actual functional 
dependence of the divergences on the (auxiliary) field background dictates the form of the necessary counterterms. The 
main issue to be demonstrated is the mutual consistency of the different counterterms required to cancel the divergences 
present in different equations. This can be done the most conveniently with reference to a unique counterterm functional. 
Also the elimination of different subsets of the variables can be realised most efhciently by manipulating the expression 
of the effective renormalised action rather than on the level of the field- and propagator-equations. 

In Section 2 we shortly outline the construction of the effective action functional starting from the Dyson-Schwinger 
equations, by explicitly formulating the closure condition. In Section 3 the leading order renormalization will be 
presented, which illustrates through a "textbook" example our logic to be followed at NLO. The main results of our 
paper are contained in Section 4, where all pieces of the NLO counterterm functional are obtained. Here, we rely 
on some detailed considerations presented in [APPENDIX Al In Section 5 we collect the pieces of the counterterm 
functional into a unique expression. The elimination of the auxiliary field in Section 6 leads to the renormalized NLO 
functional written exclusively in terms of the original fields and their propagators. Alternative elimination of the 
propagators of the pion and sigma fields provides us with the NLO renormalized effective potential as function of the 
vacuum expectation values of sigma and auxiliary fields. Some details of this procedure are given in [APPENDIX Bl 
We summarize the extensive material of the paper in Section 7. 



II. FROM DYSON-SCHWINGER EQUATIONS TO THE GENERALISED EFFECTIVE ACTION 

At the level of the generating functional Z[J] one introduces the auxiliary field a into the 0{N) model through the 
functional Hubbard-Stratonovich transformation 13011: 



-a'^(x) + -\ —a(x)(f)'^(x) 
2 ^ 2 V 3A^ ^ '^^ ^ 



cx exp <idx 



24N 



where cjp — + n'^. Then, the extended action in constant background ^/Nv reads: 



S[a, 7r„,a,i;] 



d-'x 



lid.aix))' + i(a^7r„(x))' - '^[a'{x)+7rl{x) + 2^va{x)+Nv 
—a^x) + iyAa(x)(a2(a;) + 7rl{x) + 2VNva{x) + Nv 



(1) 
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The infinite hierarchy of Dyson-Schwinger equations is obtained from the master equation [31[ 



Ja{x) 



5(I)a{x) 5(I)a{x) 



X 



(2) 



by functional derivation with respect to the field. The notation on the right hand side means that in the derivative of 
the action a generic field 4>b{x) is replaced by (j)B{x) — 4'Bix) + Gbc{x, y) s,pc(y) ^'^'^ resulting expression acts on 
the unity in the field space. Ja{x) is a current coupled to the field (j>A{x) while r[(/)B] is the effective action from which 
the inverse of the exact two-point function of two fields (pA and (f>B can be derived as iG^g{x,y) = S'^T /S(j)A{x)S(j)B{y)- 
The index A denotes the "flavor" of the field {a, a or 7r„) and possible 0{N) indices. Repeated capital indices mean 
summation over the range of "flavor" indices. 

Since Ja is the current for which the field expectation value of has a prescribed value, if one shifts the auxiliary 
field a{x) by its rescaled "spontaneous" expectation value, that is a{x) — > a{x) + a^/SN/X, then in ^ one has to take 
Ja — > for A e {a, a}. One obtains in this way the quantum equation of state for v and the saddle point equation 
for a : 



Nv \—rr? + id] + i^J^Gaa{x, x) — 0, 

— r " + i [^^' + - x) + Gaa{x, x)] = 0. 

A Z 



(3) 



Here, Gj^jGaa and Gaa are components of the exact propagator matrix defined in the "flavor" space spanned by the 
fields a,cr, 7r„. (For the sake of brevity double indices are used only for the mixing components of the propagator 
matrix) . 

The next layer of the Dyson-Schwinger equations can be written for the propagators by differentiating ^ with 
respect to the field and setting in the end Ja — > 0. One obtains: 



i{G^^)aa{x,y) 
i(G'^)aa{x,y) 



3N 

1 a 

2 V 37V 
3'^ 



GaA{x, z)GaB{x, w)VABa{z, W , y), 
GaA{x, z)GaB [X, w)TABa{z, W, y) , 

Gc,A{x,z)GaB(.X,w) + Gt,^a{x,z)Gt,^b{x,w) V ABa{z,W,y), 



iiD^^)aaix,y) 
i{D^^)aaix,y) 

i5nrnD~'^{x,y) - \j ^ I I G aA{x , z)G.^,^B {x , w)T AB7r^{z , W , y) . 



The tree-level propagators appearing here have in Fourier space the following expressions: 



where we introduced the shorthand notation 



-1, 



iiD ^)acrik) = i\l -V, 



(4) 



(5) 



Exact 3-point vertices denoted by T ABc{x,y, z) = 5^T / 54>A{x)&4>B{y)54>c(.z) also appear in these equations. The 
infinite set of DS-equations can be closed in a simple way still treating the one- and two-point functions dynamically 
by setting for these vertex functions their tree-level (classical) expressions: 



(6) 



Taaa{x,y,z) = i\j -^^{x - y)S{x - z), r^^^^a{x,y,z) = SnmTcrcra{x,y, z) . 



For any other set of indices the 3-point vertex vanishes at classical level. 

With this closure of the Dyson-Schwinger hierarchy, the set of equations given in ^ and (H]) can be derived upon 
variation with respect to the corresponding quantities from the following multivariable generalised effective action: 



r[a,v,GT,,g] 



(m^ - ia) Nv^ + 



37V 
2A' 



{N - 1) {\nG-\k) + D-\k)G^{k)) + Trlng-i(fc) + Tr {V-\k)g{k)) 
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+ i 



A 



UN 



G„„(fc)((iV - l)G^{p)GAp + k) + GUp)C{p + k)) 
Ar[a,v,G^,g]. 



' k J p 

+2Gaa {p)Gaa {k)Gaa {p + k) 



(7) 



Here, and A are the renormalized couplings, and Ar[d, w, Gjr, ^J] is the yet undetermined counterterm functional. 
Q and V are two symmetric 2x2 matrices with components Go-o-jGo-q, Gaa and D^^,D„ai D^a, respectively. Their 
inverse matrices are denoted by and 'D~^, respectively. 

The functional in ([7]) corresponds to a specific ^-derivable action. Except for the last term in the square bracket of 
the last integral, it reproduces pieces of the NLO 2PI effective action presented in (44) and (55) of [32]. This additional 
term, which comes from the self-energy of the second equation in is catalogued as next-to-NLO (NNLO) in [33 | 
and corresponds to Fig. 7a of that reference. 

Having in mind the 0{N'^) determination of the effective action we realize that most terms of the square bracket 
of the last integral in ([7|) contribute only at NNLO. Therefore we truncate further the functional in ([7]) and will work 
with the approximate effective action 

T[a,v,G^,g] = i (m^ - id) iVu^ + -^a^ 
z 2A 



+ i 



12 



(N - 1) (InG-i(fc) + D-\k)G^ik)) + Trlng-\k) + Tr (p-^ (fc)g(fc)) 
/ / Gc.^{k)G^{p)G^{p + k) + AT[a,v,G^,g]. 



(8) 



This truncated form of the functional influences the propagator equations of the a — a sector (at LO) and of the pion 
fields (at NLO). We emphasize, that even if one preserved the omitted terms, the 0{1/N) solution of that approximation 
would not correspond to the full 1/iV-accurate solution of the Dyson-Schwinger equations in the a — a sector. This 
is a consequence of the closure ([B]), which does not take into account that the vertex Taaa starts contributing at 
order through a one-loop diagram in which the pion multiplicity compensates the suppression resulting from 

having three TT^aa vertices. This results in contributions of 0{1/N) to both i{G~^)a-a and i{G^^)aa through the terms 
GaaGaa^aaa and GcraGa-oJ^aaa of the sccoud and third equations of (HI). 

Renormalizability of this approximation will be investigated at the level of the derivatives of ^ up to the next- 
to-leading order of the large-TV expansion of the equations for the pion propagator and the background fields. One 
attempts the construction of an appropriate counterterm functional Ar[d, v, G^r, g]- This countertem functional allows 
for a uniform treatement of the counterterms and also makes possible to keep track of the effect a counterterm determined 
from one equation has in the renormalization of another equation. Some new insight will be offered when compared 
to other approaches, where the propagators are not considered as variational variables 20, 2lj, or when the quantities 
related to the auxiliary field are eliminated [23|. We shall demonstrate the validity of Goldstone's theorem also for the 
renormalized NLO propagators. We start our program with the short description of the leading order analysis. 



III. LEADING ORDER (LO) CONSTRUCTION OF THE COUNTERTERMS 

A. Saddle point equation (SPE) 



Taking the derivative of the functional in 



with respect to a one arrives at the expression: 



Sa 



[a,v,GT,,g] 



3N 



N 



+ / D^{k) +c.t 



(9) 



where we replaced Gjr, originally appearing in the integral above, by Dj^ introduced in The last term denoted by 
"c.t." is the contribution of the counterterm functional which has to be constructed to ensure the finiteness of this 
equation. The structure of divergences for the tadpole integral above is given in lAPPENDIX Al Using one can 
see that the finiteness of (|9]) for any value of a is ensured by the following 0{N) counterterm functional: 



Ar 



N 



T, 



(2) 



Mo 



4" ' ^• 



(0) 



(10) 
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Here, ^ and rj°^ are the quadratic and logarithmic divergences of the pion tadpole as given in (|5D|) and ([5^ of 
lAPPENDIX Al Then, from ^ one obtains the finite saddle point equation 



3iV N 



a-i— [v^ + T^) =0, (11) 



where is the finite part of the pion tadpole integral. 



B. Equation of state and Goldstone's theorem 



The leading order term in the derivative of the functional F with respect to v is of order N 

^[a,v,G^,g]=NvM^. (12) 

The right hand side is finite in itself, it does not necessitate the introduction of any 0{N) counterterm. 

Since the equation of state is obtained by equating the r.h.s. of p2)) to zero and the leading order inverse pion 
propagator is of the form ([5|), one immediately sees that Goldstone's theorem {D~^{k = 0) = 0) is obeyed. 



C. Leading order propagator matrix in the a — cr-sector 

The only entry of the 2x2 inverse propagator matrix which receives LO correction is (G~^)aa ■ 

i{G-')Uk) = -1 + ^I^{k) + ^Tf + c.t. . (13) 



The definitions of and of its finite part in terms of D.^ , which at LO replaces , are given in lAPPENDIX Al 
(see (|68p ). In order to make this equation finite one has to introduce another piece into the counterterm functional: 

Ar"" = ^Tf |^G„„(fc). (14) 

The LO matrix elements of the 2x2 propagator matrix are then 

GS"j(fc) = (^1 - ^4^(fc)) G(fc), G(°2(fc) = -^D-\k)G{k), G^^l{k) = ^^|\vG{k), 



(15) 



where 



^G-'{k)^{k^-M') (^l-^4^(fc))-^.2_ 



(16) 



In the broken symmetry phase all elements of the LO propagator matrix have common pole structure determined 
by G(/c). This is the manifestation of the hybridization for w 7^ of the longitudinal field component a and of the 
composite field a ^ a'^ -\- tt'^ [ssj . This feature is relevant when studying dynamical aspects of the phase transition at 
finite temperature. 

Using the first entry of (|15p one can derive the following relation between the LO sigma and pion propagators: 

G<".W.Z,.,.)-4„^5ffig^. (IT) 

This relation will be very useful for the divergence analysis at NLO. 



IV. NEXT-TO-LEADING ORDER (NLO) CONSTRUCTION OF AF 

The construction of the NLO counterterm starts with discussing the l/N expansion of the pion propagator. Its 
asymptotics will be completely determined by certain integrals of the LO propagators. We shall see that the counterterm 
functional can be determined by the asymptotics of the LO propagators and of the NLO self-energy of the pions. In 
order to explicitly demonstrate the NLO renormalizability for all values of v and a we investigate the derivatives of the 
functional ([5]) with respect to these variables. The mutual consistency of the counterterms renormalizing these three 
equations is fundamental for the outcome of our analysis. 
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A. Pion propagator, equation of state and Goldstone's theorem 



The pion propagator at NLO in the 1/N expansion is given by 

^G~\k) = ^D-\k)-^^ J^Gl°l{p)DAp + k) + c.t. . (18) 

Since we need the pion self-energy to 0{1/N) accuracy, we replaced Gaa with Gal and G^r with D^r in the above 
integral. 

In order to determine the counterterm contribution in (jlSp one has to study the divergence of the NLO self-energy. 
As shown in [APPENDIX Al this divergence is momentum independent, i.e. there is no need for infinite wave function 
renormalization, in accordance with [l2|. Using the first two entries of (jlSp together with (|17p one can write 

G(o) („) = I ^ ^""(P) (19) 

Then, one sees that the momentum independent divergence is determined by the first term of (|19p : 



div 



pl-\lP{p)/6 



Tdw{M'). (20) 



div 



This divergence is worked out explicitly in [APPENDIX Al with the result 



fdiv(M2) = ri^) ~ ^(Af2 - M^)Ti^\ (21) 

where the quadratically and logarithmically divergent integrals ri^' and Ti^-* are defined in ([71)1 . 

The second term of (fT9)) does not contribute to the divergence of the first integral of (pO)) . since iterating (|17p once 
one recognizes that by logarithmic(!) power counting the following integral is actually convergent: 

7. (22) 

{l-XI^{p)/6f 

Due to this fact we do not encounter any divergence proportional to v'^ in the NLO pion self-energy. 

It is instructive to point out here a peculiarity of the resummation procedure as compared to the order-by-order 
perturbative renormalization. Namely, when in the above integral the denominator is expanded in powers of A then at 
nth order of the expansion one finds a A"(log A)""*"^ type divergence. Through formal resummation of this divergent 
series a finite result is obtained. This argument explicitly shows that in a resummed perturbation theory the structure 
of the counterterms can be different from that seen at any given order of the perturbation theory. The same effect was 
noticed in fl3| in connection with the wave function renormalization constants of pion and sigma fields which arise from 
imposing renormalization conditions on the residua of their propagators. At NLO in the expansion they are finite 
whereas in an expansion to any given order in the coupling they appear divergent. Similar consideration applies to the 
divergence appearing in (j20p . Expanding the denominator of the middle expression into powers of A one finds the usual 
quadratic and logarithmic divergences characteristic for the perturbative contributions. Resummation modifies these 
divergences as can be explicitly seen in the last two formulas of [APPENDIX Al 

Since the necessary counterterm in p8|) compensates Tdiv, one readily finds the required counterterm functional piece 
upon functional integration of (pij) with respect to G-^ and multiplying it by A/6: 



AP- = -- 
6 



r(2)-^(M2-Mo2)Ty) 



G.(fc). (23) 



fe 



Next, one investigates the renormalization of the derivative of the 2PI effective action with respect to the background 
V. At NLO in the 1/iV expansion this is given by 



^ Nv M^ + — G{k) + c.t. , (24) 

where for the second equality we used the last entry of (|15p . The counterterm functional AP'^ determined above does 
not contribute since its derivative with respect to v is zero. 

The equation of state is obtained by equating the r.h.s. of (j24p to zero. Its unrenormalized expression obviously 
implies when confronted with p8)) the validity of Goldstone's theorem with 0{1/N) accuracy. Note that, as it is 
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well known, Goldstone's theorem is not followed if one proceeds in strict 2PI sense which requires the self-consistent 
determination of the full propagator without expansion in 1 /N. 

In a renormalization procedure which preserves Goldstone's theorem one should construct a counterterm which 
does not depend on Gtt, therefore does not interfere with its already renormalized equation. Since the divergence in 
([M)) coincides with the divergence of the NLO pion self-energy, the new contribution to the counterterm functional is 
obtained upon integrating with respect to v the expression given in (|2ip multiplied by Aw/3. One finds 



6 



(25) 



We conclude this part by giving the finite pion propagator at NLO in the expansion, including also the contri- 
bution of the counter- functional AF'^. It reads as 



37V 



(fc) = ^ / G(^l{p)DAk+p) - TMM'). 



(26) 



A remarka.ble feature of this 0{1/N) renormalized solution is that it satisfies Goldstone's theorem for arbitrary values 
of a ! 



B. Saddle point equation 



In writing down the derivative of the effective potential with respect to a one has not to forget about the contributions 
of the (5-dependent counterterms AF"'^, AF'^, and AF" constructed above (see (|10p . (1^51) . and 



da 



N 



.N 



T. 



(2) 



(M2 



G,(fc) 



12 " 



(C(fc)-G,(fc) 
^ ' ^ G^k) 



(5AF"^0( 



a 



5a 



(27) 



The contribution of the counterterms determined by the renormalization of the equation of the inverse pion propagator 
(US]) and of the equation of state is displayed in the last but one term. 

The last, yet undetermined part of the NLO counterterm functional, i.e. AF"'", provides the NLO completion to 
AF"'^. In order to determine it one first has to evaluate with NLO accuracy the pion tadpole in the second term of 
the r.h.s. of ([77]) . Taking the inverse of G^ given in (pS)) and expanding it to 0{1/N) one obtains 



G.(fc) 



D^k) 
DAk)-i 



A 

3iV 
. A 

37V 



Dl{k) / G(S(p)I?.(fc+p) 



J(m2) - fdiv(M2) / Dlik) 



37V 



iTdiv(Af') 



Dl{k) 



where for the second equality we used (fT9|) and introduced the following functions 



J(7V/2 



Dl{k) 



D^p + k) 



J{M 



/r2\ 



Dlik) 



DAp + k)G^^}{p) 

; {i-xi^{p)/&f 



lpl-XI^{p)/6' 

Collecting all C'(7V'') (NLO) divergent terms in ((27)) one sees that SAT°'''^{a)/5a is determined by 
(5AF"'0(q,) 



(28) 



(29) 



Sa 



Jdiv(M2)-fdiv(M') 

G^°}{k)-DAk) 



Diik) 
A2 



A^ 

i—v'^JdUM^ 

io 



+ 12 



D^k) 



(30) 



where jdiv(Af^) and Jdiv(7Vf^) denote the divergences of the integrals defined in Note that to the order of interest 

it was again allowed to replace in the last two terms the original G^ and G^a by Dt, and G^aa , respectively. 

The important question of consistency inquires whether the previously constructed counterterms cancel all subdi- 
vergences of J{M^) and the w^-dependent divergence of second and third terms in ([30]) . The double integral J{IvP) 
has only overall divergence, since both k and p integrals are individually finite. This divergence is determined in (|78p 
of I APPENDIX Al The divergence of the third term of ((30)) is determined in ((76)) and also shown to be proportional 
to . One finds that the sum of these two w^-dependent divergences is canceled by the u^-dependent counterterm 
contribution appearing in the last term of ((50)1 . 

The detailed T = analysis of the divergence structure of Jdiv(7V^^) given in [APPENDIX A) is based upon the 



FIG. 1: The appearance of vertex type subdivergences in the first double integral of (I29p as illustrated at leading A order in 
the expansion of the denominator of the integrand. The factor of three indicates the possible ways in which two lines of the 
setting-sun diagram form the bubble which corrects the vertex at 2-loop level. The cross denotes the associated lowest order 
counterterm. 

explicit expressions of some integrals. The presence of subdivergences is reflected by the appearance of divergent 
terms proportional to ln(M^/MQ). These should cancel if the approximation is renormalizable, since in this case only 
divergences proportional to powers of a (that is M^) are allowed. The result given in ([80| for Jdiv(M^) shows when 
combined with the second term of the square bracket of ((30|) the cancellation of the subdivergence Tdw{M'^)lTr{p — 0). 
The cancellation of this self-energy type subdivergence of the double integral J{M'^) is expected in view of (pS)) . 

The double integral Jdiv(Af^) has also a vertex type subdivergence as illustrated in Fig. [T] at leading order of the 
expansion in A. This is canceled, as it should, by the last integral of (j30p . One can see this analytically by separating 
in the difference of the two terms of the square bracket of (jSOp a contribution proportional to T^^T^ which on its turn 
cancels against the contribution of the last tadpole integral in ([30| . 

With all subdivergences and w^-dependent divergences of ([50]) canceled, SAr°'''^{a)/6a is determined by the overall 
divergence of J{M'^) and that of the last tadpole integral. Its expression reads: 



Since the above expression depends only on a it will have no "back-reaction" neither on the propagator equations nor 
on the derivative of the effective potential with respect to the background. The appearance of terms proportional to 
l/(167r^) reflects some allowed arbitrariness of the subtraction procedure. 

We do not present here the renormalisation of the partial NLO corrections which would occur in the a — cr sector 
should one use the effective action functional ([7]). All counterterms necessary for the renormalisation of these pieces 
are of 0{\/N), not contributing to the renormalised effective action at NLO in the 1/A'' expansion. 

In conclusion, only divergences proportional to zeroth or first powers of remained which upon integration over a 
determine the d-dependent counterterm functional /!^T°''^{a). The counterterms induced by the renormalization of the 
NLO pion propagator and of 6T /5v played an essential role in the cancellation of subdivergences and of w^-dependent 
divergences present in the exp ression of 5T /5a. No limitations whatsoever showed up on the value of and/or w^, in 
contrast to the findings in [lo, HH ■ We shall return to the discussion of this difference after reducing the effective action 
to the effective potential depending only on the background fields in [APPENDIX Bl The functional corresponding 
to the NLO l/A^-expansion of the classical vertex approximation (BVA) is renormalized without any constraint. Our 
result shows that this can be achieved also without introducing unconventional counterterms . 



In this section we collect into a unique expression AP all individual pieces determined in Eqs. PH)) . pi]) . (^5]) . (pS]) . 
pip , and express it in a conventional form, in which one associates them with the renormalization of different couplings 
appearing in the terms of P[d, v, G^r, Q] (Eq- (©)• The counterterm functional reads: 



5AP"'0(a) 
So. 




A^ - Af2" 
12 



(31) 



V. 



THE EXPLICIT FORM OF THE COUNTERTERM FUNCTIONAL 




(32) 



where the counter-couplings are given by the following expressions: 
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— — 



■(3m2 - M^) 



487r2 



(33) 



It is interesting to note that the term proportional to Ski in (|32p has no correspondent in ([5]). Moreover, the terms 
proportional to Skq and 6k2 correspond in ([5]) to terms not proportional to any coupling of the original theory. Rather 
they renormalize numerical coefficients related to the way the auxiliary field a is introduced by the parametrization of 
the Hubbard-Stratonovich transformation (cf. ([T])). They correspond to different possibilities of defining the two-point 
function of the auxiliary field. They are calculated to different orders in 1/N since the corresponding terms contribute 
at different levels of the expansion. After scaling back the fluctuating part of the field a to a the two counter couplings 
do agree to leading order in 1/N. This feature is reminiscent of the renormalisation of the operators corresponding to 
different definitions of the n-point functions in 2PI-approximations. The terms with dg can be regarded as renormalizing 
the coupling g = ^/X/(3N) through which the auxiliary field couples to fields of the 0{N) multiplet. One interprets 
the renormalisation of all these operators not occuring in the original model as renormalisations of the two variants of 
the 2-point functions of the composite field 0^ and of the ip'^aa vertex (see ch. 30 of (33|V 

It is notable that Sk2 determines the renormalization of the coupling A following the formula 



_6_ 



6 4, 6 
A + iv'^"^= A 



N + 8 
N 



^(0) o—T^i) 



6 647r2 



(34) 



This can be seen by looking at the terms proportional to in the classical part of the functional. This relation is only 
intermediary, since it will change in the course of the elimination of a. The part of the counterterm Sk2 proportional 
solely to r^^"* appears already in the literature and follows the one-loop /3-function of the model [sHi • The entire NLO 

part of Sk2 proportional to Ta^'^ is missing in the analysis of [13, [lH (see e.g. Eq. (23) of [21.] for the expression of their 
counterterm bi). 

Introducing the following notations 



Sk2 = N5k^°'^ + 5k. 



.(1) 



As = A + 5X^, 5X^ = SX^°'> + -5A(,i) 



N 



one readily obtains 



T, 



(0) 



xf = x + sx(^^ 



X 



I + XTf^/Q 



SXi' 



(o)^ 

B ) 



'Sk, 



(1) 



(35) 



(36) 



Comparing with Eq. (25) of [2^ one observes that 6Xa^ is the counter-coupling of the 0{N) model formulated with 

its original variables and considered at LO in the large- iV expansion. Likewise A^^ is the bare coupling of the model 
in the large- iV limit. We shall see in the next section, where the auxiliary field will be eliminated, that at NLO in the 
large- iV expansion the bare coupling of the 0{N) model differs from As, it turns out to be a combination of As and 
6g (see ge]) and (|47)) ). 

With the counterterm functional AF explicitly determined by (j32p , one can give now in a compact form the functional 
introduced in 

rNLo[a,v, G^,g] = {ml~icay+iSKia+—a^-^ [Tr In G - \k) + {N - 1) In G-\k)] 
2 2Xb 2 J I, 



lj^[e-ml+ ^ca] [{N - f )G.(fc) + G(°j(fc)] - ^v^|\ G^^l{k) 



where As is defined in (|34p and we introduced the following notations: 



m% = + —6m'^, c=l + —Sg, c = i5ko, Ti{k) = -i j G^{p)G-„{k + p). 



(37) 



(38) 



VI. ELIMINATION OF THE AUXILIARY FIELD 



In this section the 0{N^) accurate renormalized functional will be established for the original formulation of the 
0{N) model by eliminating the auxiliary field a and the propagator components related to it (i.e. Gqo- and Gaa)- 

In order to achieve this one substitutes into (|37p the expressions of a,G'^aa, and G^aa as found from their respective 
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equations. 



A. Determination of rNLo[i', G,r, Ga] 



For rewriting the terms depending on Gal and Gal one exploits their representations which allow the expression of 
the result fully in terms of n(fc) and ■ The latter will be replaced with Go-, the exact sigma propagator of the 0{N) 
model. In this way one finds 



-^Trln(?-i(fc) = -\ j\ndoig-\k) = -'- j\n 
In going from the first to the second line above we used that to 0{\/N) accuracy 



ln(zc) + 0( - 



I / G^(p)G^(k+p) 



\ I D^{p)D^{k+p) 

J V 



T. 



(0) 



d ' 



(39) 



(40) 



and therefore one has 



^i^{k) = c-^h{k)^c-^n{k) + o 



where the neglected 0{1/N) contribution is finite. 

Using (fT9|) for G^aa one writes the following chain of expressions for the last term of ([37]): 
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jc^^likmk) = -| jG^Sl{k) + \ jG^^l{k){c^-^U{k) 



6c 



G.(fc) 



/ i + o(i 



N 



(41) 



(42) 



In writing the second line above one again uses (|4T|) . The first term on the r.h.s. above is canceled against the last but 
one term of (|37p . Finally, for the last term in the second line of (|37|) one uses (|4T|) to write: 



(43) 



As a short digression from our current task we mention that one could proceed by further eliminating also the pion 
and sigma propagators using their respective NLO and LO equations. Then one obtains the renormalized version of the 
effective potential as function of a and v studied in [l^, [2l| . A sketch of this derivation is presented in lAPPENDIX Bl 
together with a check of the renormalization of the saddle point equation coming from this potential. 

Now, we proceed instead with the elimination of a from (j37p keeping the variables of the original theory, i.e. w, G^, 
and Gct- Actually, the simplest way is to complete to full square the functional depending quadratically on a and then 
to use the saddle point equation 5T/6a. — 0. Replacing G^ct by the exact propagator Ga of the 0{N) model one obtains 
from the d-dependent part 



24iV 



Ncv"^ +c [(iV - l)G^(fc) + Ga{k)] - 2Ski 



1 2 



(44) 



Putting together all above pieces one also uses that in view of ([36]) A/c = A^-* and obtains 



N 

r[?;,G^,G^] = y I "I 



2 XbcSki 
^ 3iV 

2 I 



N 



Xbc' 

'I 

24 



[(A^-l)lnG;i(fc)+lnG;i(fc)] 



m 



XbcSki A^c^ 2 



B 



3N 



[(iV-l)G,(fc) + G,(A:)] 



^„2 

6 



+:^(iv-2)( / GAk)] + 



GAk) 



k 1 - A^°^n(fc)/6 



\(0) 

In ( 1 - -^n(fc) 



"T2~ 



G.(fc) / G,{p), 



(45) 
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where we omitted terms of order 0{1/N) and a divergent constant ^ Snf coming from (M)) as well as the irrelevant 
divergent last but one terms of ([39)l and (|42|) . 

The bare couplings appearing in (|45p have to be given with an accuracy corresponding to that of the functional ([8|) , 
therefore one writes the counterterms as a sum of the LO and NLO contributions: 



X + SX, 



<5A = a(") + l(5AW, 



2 XbcSki 
rriB 



2,N 



6m 



2(0) 



1 . 2(1) 

— dm 

N 



(46) 



With help of ([M]) . ([55)) . ([55]) and the separation (5ki = N6k^^^ +(5k^^\ the counter-couphngs above can be given 
in terms of the counter-couplings of the model in the auxiliary field formalism as 



<5A(i) =(5AW+2A^"'<5.g, 



Sm 



2(0) 



(O)x^(o) 



-X'^'Sk- 



X 2(1) X 2 1 

dm ~ dm 

3 



SXi'^S^^°^+xf [5n['^+S^^°^5g 



(47) 



Using ([^5]) in ([IS]) one obtains at NLO in the 1/iV expansion the renormalizcd functional of the original 0{N) model, 
that is without the auxiliary field, in the following form: 

Tnlo [v,G^,G,] = -'-j^[{N-l) (In G-\k) + &-\k)G^ (k)) + In G'^k) + &-\k)G, [k)] 



iV , , X . X 
2 24 24 



A 



(0) 



GAk) 



X 



(0) 



GAk) 
Ga{k) 



X 
12 



G^(fc) / G,{p) 



k 1 - X^°'>U{k)/6 



i(0) 



^ G.(fc) In l-^n(fc) 



12 



X 



(0) 



SX 
12 



/ [(iV-l)G,(fc) + 3G,(fc)] + 



6 

[(iV-l)G,(fc) + G,(fc)] 
N 



G.{k) 



G^k) / G,{p) 



(48) 



where we introduced the usual tree-level propagators for the sigma and pion fields as 



A 



A 



i&;;\k) ^ k^ - m^ - '-v\ i&-\k)^k^ -m"- -'-^v^. (49) 
If one forgets about the counterterms, the expression in ([48]) coincides with the 2PI effective potential obtained in 



[27l |. The terms above can be combined in a way which makes explicit that in the large- expansion there are only 
two bare couplings, namely m^ + 5m^ and A -I- 5X. Would one attempt a selfconsistent solution of the equations arising 
from the variation of Fnlo this would not be true (see the concluding part of this subsection). We emphasize also that 
since the above functional is 0{N^) accurate, in terms involving the sigma propagator G^ one does not need the NLO 

part of the counterterms, i.e. (5m^' ' and 5X^^^ . When differentiating with respect to G^r one should remember that also 
n(A:) is a functional of G^! 

The interpretation of the terms in (j48[) , which makes explicit the infinite series of diagrams summed up in the present 
treatment is as follows. The first two and the last two lines represent the 2PI effective potential of the 0{N) model 
at Hartree level of the truncation and at NLO in its large- A'^ expansion. The remaining four terms represent the NLO 
contribution of the 2PI vacuum diagrams beyond Hartree level. The u^-dependent part of these terms can be rewritten 
as 



A + 

6 



G.(fc) 



(A + aW)n(A:)/6 _A + a(o) 



1- (A + (5A(o))n(A:)/6 



' / G.(fc)^ 



A + (5A(") 
6 



Uik) 



(50) 



The ii^-independent part can also be written as a sum: 



A + JA(") 
6 



n(fc) 



In ( 1 U{k) 



E; 

n=2 



1 fX + 5^ 



U{k) 



(51) 



It is easy to show that not considering counterterms, these terms correspond to the two sets of diagrams given in 



12 



FIG. 2: The two sets of vacuum diagrams which contribute beyond Hartree level and at NLO to the 2PI functional. Solid 
(dashed) line represents sigma (pion) propagator, while wiggly line represents the background i;. The dots indicate any number 
of pion bubbles. 




FIG. 3: Diagrams corresponding to the n — 2 term in the sum on the r.h.s. of (|50|l . The cross represents the counter-coupling 



Fig. [21 Counterterm diagrams corresponding to the n — 2 term of the sum in (|50p are displayed in Fig. [3) Similar 
diagrams with different number of pion bubbles can be drawn for the other terms appearing in the sums in (|50p and 
([5T|) . A direct way to obtain consists of summing up all these diagrams with the associated combinatorial factors 
determined by the Feynman rules. 

An interesting question is how the counterterms of the last two lines of (^S)) are related to the set of counterterms 
introduced in [ij]. That structure was fixed by the possible 0{N) invariants of the model. Comparing the expressions 
presented in Eq. (6) of [1^ with those of the last two lines of one obtains the following unique relation for all of 
them: 

SX4 = SX^ = SX^ = (5Af = 6X^ = 6X, 6ml = ^^^l = Sm^- (52) 

This result is actually expected on general grounds, since in the strict expansion used in our approach different 
definitions of the two- and four-point functions coincide. The above coincidence represents rather a check of the 
correctness of our NLO renormalization. 



B. Renormalizability checks on rNLo[i', G,r, Ga] 



The significance of (pg)) is that it displays all counterterms which guarantee the renormalizability of the resummation 
of the perturbative series, a resummation induced by the large- expansion. The finiteness of the equation of state 
and the self-energies obtained from its respective variations is ensured automatically, since this feature is "inherited" 
from the finiteness of the same quantities achieved in the formulation with the auxiliary field. Still, it is an instructive 
exercise to check this feature directly. Exploiting the structure of our previous analysis done in the auxiliary formulation 
of the model we shall show the finiteness of the equations directly obtained from 

The inverse pion propagator at NLO in the 1/N expansion is given then by 



where the nonlocal and local parts of the self-energy are 



(53) 



Xv' 



1-XIIf{p)/6 3 (l-AnF(p)/6) 

Ar / O /" „ ^ \ Xb 



G,(fc+p)-Tdiv(M2), 



GAk) 



6N 



Since the local part has LO and NLO contributions one writes 



3A^ divV j 



For the nonlocal part we used that B/A^'' — IT = 6/A — Tip/Q and Tjiv is given in ([H 



(54) 



^2(0) _^ M"^^^^ /N and expands the pion 
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propagator to 0{\/N). With help of the integrals defined in (|29p one obtains 



M 




X 



(0) 



Si 



1-^^^(0)1- — 



6 



A 



v' + I D^{k) ) , 
fc 



r 2(1) , ^rf, (T,r2'fi)\ 

dm ' ' ■■ ' " ' ' 



J(Af2^"Vrdiv(M2(°)) / Dl{k) 



Ga{k)^D^{k) 



. A^ 2 -r/, ^2(0)x ^f^A*^ 
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2A 



(0) 



(55) 



where in comparison to its definition given in ([5]) depends now on M^*'°''. We used also that to leading order 

n(p)-/,(p). 

The equation for Af ^ in ([SS]) is the usual gap equation at Hartree level of truncation of the effective action. This 
was analyzed in [2^ and the counterterms which can be determine from this are 5\^'^^ and Jm^'"' given in (|47p . 

Observing that the left hand side of the equation for Af^^^'' in (|55p is finite one obtains the following relation between 
counterterms and divergences: 



5m 



2(1) 



^Tdiv(Af 



2(0)n 



Jdiv(M 
i(5A(i) 



2(0)x 



fdiv(M^(''^) 



Dl{k) 



A^ 



2(0) 



2A 



(0) 



D^ik)] + '- {G„{k) ~ D^{k)) 



(56) 



div 



Using the divergence analysis of I APPENDIX Al one has all divergences and integrals expressed in terms of Af^^'^'' for 
which one can substitute its finite gap equation K'p''^^ — m^ + A/6(ti^ + T^). Requiring the vanishing of the coefficient 
of + determines 5\^^\ while the remaining overall divergence determines 5m^''^\ Both are in accordance with 
(gZl). 

The equation for the inverse sigma propagator obtained from (|48p is 



iG-\k) 



Af2(°) 



k(O) 



Ag^/.(fc)/6' 



(57) 



which is finite, since 

A^2(0) is finite and 6/ A 



v(0) 



6/A-/f/6. 



Using the relation between the LO pion and sigma propagators one can show that the derivative of 
to V reads as 



5T 



NLO 



5v 



[v,G^,G,]^^NviG-\k = Q), 



with respect 
(58) 



which is also finite, since we showed that G~^{k) is finite. It also displays the validity of Goldstone's theorem. 

We close this part by mentioning that the LO sigma propagator equation, the NLO pion propagator equation and the 
equation of state derived from (j48p can be obtained also with the Dyson-Schwinger formalism of Section 2. Concerning 
these equations, the large- TV expansion closes the hierarchy of DS equations at the level of complete LO renormalized 
r^^o- and r^^TTTT vertex functions, which include also one-loop pion corrections. Some details can be found in [36| for 
Gu and F^^o- (see also [33| for the relation between the truncation of the Dyson-Schwinger equations and the 1/7V 
expansion). This means that our investigation implies also the renormalizability of the Dyson-Schwinger equations at 
NLO in the large- expansion. 



VII. CONCLUSIONS 



We studied the renormalizability of the 0{N) model at next-to-leading order in the large- A^ expansion, at zero 
temperature. We constructed the 0{N^) counterterm functional of the model in auxiliary field formalism by studying 
the renormalization of the derivatives of a generalized effective action functional with respect to its variables. Expanding 
the propagator equation of the pion to 0{1/N) in the large- A^ expansion we showed that the renormalization can be 
achieved for arbitrary values of the background and auxiliary field, in a way that respects the internal symmetry of 
the model (e.g. Goldstone's theorem). This can be expected on theoretical grounds, since divergences are determined 
only by the asymptotic behavior of the propagators and because any consistent resummation of the perturbation 
theory should resum also the counterterm diagrams associated to the perturbative series. Although, one can anticipate 
the consistency of the auxiliary field technique and of the large- A^ expansion in dealing with perturbative series, the 
difficulty we face when trying to infer the renormalizability of the model in a given approximation from the fact that 



14 



the model is perturbatively renormalizable is that one cannot easily keep track of what partial series of counterterm 
diagrams is actually resummed at a given order of the large- A'^ expansion. In consequence, the actual analytic check of 
the renormalization of a given approximation is unavoidable. 

The elimination of the auxiliary field and the related propagators, while keeping the dynamical sigma and pion 
propagators, makes transparent the classes of diagrams containing also counterterms which are resummed in the 
original 0{N) theory at NLO of the large-iV expansion. The explicit form of the counterterms is given here for the 
first time for the theory using auxiliary field and also for the original formulation. In the original theory the action 
functional contains only two counterterms, a coupling and a mass counterterm, both having LO and NLO parts. The 
propagator equations and the equation of state derived from that coincides with the IPI Schwinger-Dyson equations 
closed at the complete LO r7r,rcr and Tt^t^^t^ vertex functions, which includes one-loop pion corrections. 

The two examples we worked out (e.g. FnloI^^t G^, G^] and rNLo[(i, v]) demonstrate that the renormalizability of the 
broadest action functional T[a,v,GmG] implies the renormalizability of the functionals arising after the elimination 
of some subset of the variables. This result obtained at T = makes us confident that the renormalization goes 
through for nonzero temperature as well and that the counterterm functional determined here will prove helpful for 
phenomenological studies in the 0{N) model. A study of the renormalization scale invariance at NLO in different 
formulations would be helpful for such applications. The method developed here can be used also to models with more 
complicated global symmetries. 
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APPENDIX A. DETAILED ANALYSIS OF THE DIVERGENCES 



Notations 

Following the spirit of Ref. [29j we will expand the propagators around appropriately chosen infrared safe auxiliary 
propagators. Since at NLO the asymptotics is determined basically by integrals of the LO propagator D7r(p) and of 
the propagator G{p) given in (fT6|) which incorporates the effect of the resummation of the chain of pion bubbles, one 
needs two auxiliary propagators. The first one is 

Goip) = (59) 
where Mq is an arbitrary mass scale. With this propagator one defines the quadratically divergent integral 

T^'^ = / Go(p), (60) 
Jp 

and the following one-loop bubble integral 

Io{p) = / Go{k)Go{k+p) = ri°) + I^{p). (61) 

Jk 

The logarithmically divergent part of the integral above is defined as 

Tf^ = / Gl{k). (62) 

Jk 

The finite part behaves asymptotically as Iq{p) ^ In — 2 — in + O Inp) , and defines together with Go{p) the 
second auxiliary propagator: 

The integrals involving only combinations of Gq{p), Ga{p) and Iq{p) will be fully included in the counterterms. With 
help of the propagator Go{p) one can separate the quadratic and logarithmic divergence of the tadpole integral defined 
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with the tree- level propagator given in ([51). At zero temperature one has: 



D.{k) = Tf^ + (M^ - AQ Tf^ + , t! = ^ {aI^ In ^ - Af ^ + . (64) 



Absence of wave function renormalization in the NLO pion propagator 

The only term which on dimensional grounds could produce a momentum-dependent divergence in the expression of 
the NLO pion propagator (fTS]) is the first one on the r.h.s. of p9)) . It gives the integral 



To study this integral one uses the following expansion 

1 _ 1 1 P + 2p-k Y 

- tD^ip + fc) - _^ _ " p2 _ + p2 _ A. p2 _ ) ' ^^^> 

In order to prove that there is no infinite wave function renormalization it is enough to look at the appearance of fc2 
in the numerator, that is at terms n = 1, 2 in the sum. Keeping only terms up to and including ©(l/p^), but throwing 
away those which vanish upon symmetrical integration (note that I^{p) depends on p^) amounts to the following 
replacement at the level of the integrand in (|65p : 



1 1 , 4(p-fc)2-fc2p2 ^^^^ 



{p + fc)2 - M2 p2 _ 2 (p2 _ ■ 

The second term above gives vanishing contribution in (|65p due to the following property which holds for any integrable 
function f{p^) upon the use of a Lorentz- invariant regulator: 

/ d'pp.p.fip') - ^ 1 d^pp'fip'). 

The first term on the r.h.s. of (p7|) gives the momentum independent divergence denoted with Tdiv in (|20p . 

Momentum independent overall divergence of the NLO pion propagator 

To find the momentum independent divergence of the NLO pion propagator one starts with ((20)) and takes into 
account that the one-loop bubble integral behaves logarithmically for large momentum. 

The divergence of the one-loop bubble integral is chosen to be T^"'' given in ([5^ . Then, one has 



M2 

-2 + In — 2 + L(p, M) 



(68) 



/x(p = 0) = -i / I?^(fc)=Tr+/^(p = 0), 4^(p = 0) = ^ln^. (69) 



Here, the momentum dependent function L{p'^,M) which determines the finite part I^{p) can be found in Eq. (8) of 
Ref. [H]. From hmL(p,M) = 2, it results that 

1 , M2 

Since Iq{p) has exactly the same form as I:^ {p), but with Af2 replaced by Afg , for large p^ one has 

,>,_,„„._LM_^,„^,„(^). 

In the asymptotic momentum region the above expression allows us to write 

1 _ 1 , A I^ip)-Iiip) ^^f 1 



l-A/^(p)/6 l-XI^ip)/6 6(l-A/o^(p)/6)2 \pnnp, 

1 , AA.f2-A.f2 liip) , ^ / 1 



1 - Xliip)/6 ■ 3 p2 _ M§ (1 - A/o^(p)/6)2 \p2 ln2p 



-TT^ ' (71) 
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where in the last Une we used Iq (p) ^ In — 2 — in + O (p ^ hip) . The neglected terms give finite contribution in 
the last integral of (|20)) . Using there l|7l|) and 

1 1 



one obtains 



(72) 



rdiv(M^) = T^2) - (A/2 - M2)z Giip) 1 - -I^ (p) 



A 



where the following divergent integrals were defined 



Ti2)= [ Ga{p), Ty^=-^ [ Gl{p)l', 

J p J p 



ip)- 



(73) 



(74) 



Using in the remaining integral of ([75)1 that J G^{p) is finite, one obtains for T^iv the expression given in (pi]) . 



Analysis of the divergences of the saddle point equation (|30|l 

We start with the differences of tadpoles involving LO sigma and pion propagators. Using p7p iteratively one finds 

Dl{k) 



A 



6 J^l-\lF{k)/Q 



In view of (j7ip one can replace I^{k) with I^^ {k) in the denominator above. Then using (j72p one obtains 



GW(fc)-i?.(fc) 
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/ G^(fc)(l--/„^(fc) 



36 ° 



(75) 



(76) 



Next, we investigate the second double integral given in l|29p. Changing the order of integration one uses P7)) and 
the following relation which can be derived from 



J^Dl{k)DAp + k) 



1 d 



to find 



JdiviM 
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t Gl{p)li{p)=Tp. 



(77) 



(78) 



To obtain the second equality above we replaced with Iq in view of ([70]) and used ([72]) . 

The first double integral given in (|29p contains an overall divergence as well as subdivergences. By changing the 
order of integration and using (j77p . one has 



Jdiv(M2) = - 
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A 7„ p2 _ 
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(79) 



div 



where we have separated a divergence independent of I^{p), which by using a relation similar to (j72p can be expressed 



as a linear combination of T^^-* and Tj*^-* . 



Since the form of the second integral in ([79]) is similar to the last one in (PO)) . the calculation follows very closely the 
determination of T^iv Using ([2T|) and (|69)) the result is 



^J,i,(M2) = zTf + (4M2-M2)zTf) 
A 



+i ^M^tP - fdiv(A/2) 1 + -T^ 
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/.(p = 0) 



(80) 
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Cut-off dependence of the divergent integrals 



.(0) 



It is instructive to evaluate explicitly the cut-ofF dependence of the divergent integrals denoted with , j_ ^ 

(2) 

and Ta . This can be also of some practical interest when one proceeds to the numerical solution of the renormalized 
equations. Going to Euclidean space with po ip'^ and using a 4d cut-off A the first two integrals defined in ((62)) and 
([60]) can be done analytically: 
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For the other two integrals defined in (|74p we limit ourselves to an asymptotic analysis and expand the integrand for 



large k. Exploiting the freedom to omit those contributions to Ta 
scheme in which 



(!) 



which are formally finite for A —^ oo, we choose the 



(47r) 
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(81) 



Here, we introduced a = A/ (967r2). One notices that for k = AIq exp(l + 487r2/A) the denominator of the integral above 
vanishes. To avoid this non-integrable singularity to occur in the range of integration, that is for A: < A, one needs 
A < Amax — Mq exp(l + 487r2/A). That means that for A 7^ the cut-off cannot be sent to infinity, there is a maximal 
value for it, which refiects the triviality of the theory. 

Performing the integral in (jSip . and obeying this restriction, one finds 
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With the same strategy one can choose 
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APPENDIX B. DERIVATION OF THE EFFECTIVE POTENTIAL V-Nho[a,v] 

One starts from ((37|) and after using ([42]) and ([43]) one finds that in view of pT]) G^ct = Go- drops out from the 
functional, which now becomes 

37V" 2 f 

T[a,v,G,] = -{ml-icay+iSnia+—a^~^iN-l) / {lnG-\k) + D-\k)G^{k)) 

Z ZAb Z Jf. 



-- I In 



j,(0) N x(0) 

(A:2 _ ^2 ^ j ^ _ ^n(A;) j - 



(82) 



Here, the last term comes from Next, one uses ([25]) . (HU, and the definitions in ([55)1 to write the inverse pion 

propagator as 

iG-Hk) = k^-ml+ica-^j:^{k), (83) 

where S^(fc) is given by the integral of (pS)) calculated with the expression Gal taken from p^ . Using this propagator 
in the first integral of ((82)) . one easily sees that when expanding it to 0{l/N) the contribution of the self-energy drops 
out and we are left with 

^NLoictTu] = —{m'^ — ica)v^+iSKia+— — cy^ — N— / \n (k^ — m% + ica) 
2 2A_B 2 Jj. 

The radiative part of this functional has exactly the same form as the effective potential given in [1, [2^ [2l| . The 
difference in the classical part corresponds to slightly different ways of introducing the auxiliary field. More important is 
that the authors of [13, |2l| restrict their counterterm functional only to terms proportional to pieces of the Lagrangian 
which are present already in the original formulation of the model. In their form of introducing the auxiliary field this 
restricts the counterterms to those proportional to a and d^. By allowing all independent counterterms to appear which 
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have dimension less than or equal to 4 in the auxiliary field formulation one might expect to have enough flexibility to 
ensure the renormalisibility in arbitrary background ^. 

In order to demonstrate that ([84| contains all the NLO counterterms, we sketch the renormalization of the SPE 
obtained by differentiating ([84| with respect to a. Using that 11, to be taken only at LO in expansion, depends 
on a through D^^ defined in ([5]), one obtains 

3A^ - -r f 2 f T- 

= —a~idKx-i—\v + / — 2 I ■-- 

A_B 2 \ J f. k'^ ~ + ica 

^\ I (G^^Kk) - DAk)) - ^ J(M^) + ^^v'Ji^n (85) 

Here, we recognized the appearance of the expression of G^^ which can be read from ([T5)) . We used the relation ([T7)) 
and for the last two terms also (jH]) and ([^5]) . 

All we have to do is to establish the connection between ([85]) and ((27)) . the latter being already renormalized. The 
last three terms of ([55]) can be found in PT)) . if in that equation one takes into account so we have to work only 
on the first three terms of ([85]) . Using the definition of the couplings, one expands them to 0{1/N) and obtains 



3A^ , 3iV r 

- — a — iOKi = —r-o; + i — 
AB A 2 



T^'^ + (M2 - M2)ri"' +iS^,['> +2S>i^^'>a. (86) 



.(i)^oA.(i). 



Since the last two terms of the second equality above coincide with 6Ar°'''^{a)/6a of ^7\ . the equivalence between 
and (l27l) is demonstrated. 
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